
*PLEASE READ: The replication for Appendix E is at the very end of this file with asteriks so your computer won't read the commands because in order to estimate the results in Appendix E you need to download Brandon Bartel's Clustergen User Written State File net from https://home.gwu.edu/~bartels/stata Then, click on the blue link called clustergen. A new window will pop up. Simply click on the link on the right side of the page that says "click here to install." This program assumes that your cluster identification variable is called "id." If your cluster variable is called, e.g., "countryid," simply issue the following command before using clustergen: gen id=countryid See the help file associated with the program for more details. To use the program, simply type "clustergen" followed by the list of level-1 variables that you want to transform. For each variable in the variable list, two new variables will be created: (1) a between-cluster version, which is the original variable name followed by "_bw"; and (2) a within-cluster version, which is the original variable name followed by "_wi."  . To do this the long way (without clustergen), one would first have to create the between-cluster operationalization of each variable one command at a time: egen x1_bw=mean(x1), by(id)egen x2_bw=mean(x2), by(id)egen x3_bw=mean(x3), by(id) etc... Generating the within-cluster transformations would entail the following commands: gen x1_wi=x1–x1_bw gen x2_wi=x2–x2_bw gen x3_wi=x3–x3_bw etc... The clustergen program greatly simplifies this task by transforming all variables in one fell swoop. 

use E:\PartyPoliticsDataset.dta
set more off
xtset stnum year, yearly



*Interacting Macroeconomic Outcomes and Minimum Wage Gap with Party Control
gen rminwag = d.minwag2adj
gen lrminwag = l.rminwag
gen inflation1 = inflation * parcont
gen inflat = inflation
gen linflat = l.inflation
gen dinflat = inflation * demcont
gen rinflat = inflation * repcont
gen unemp1 = unemp * parcont
gen dunemp = unemp * demcont
gen runemp = unemp * repcont
gen aveincadj1 = aveincadj * parcont
gen aveincadj2 = aveincadj * demcont
gen aveincadj3 = aveincadj * repcont
gen mingap1 = d.newmingap * parcont
gen dmingap = d.minwaggapadj * demcont
gen rmingap = d.minwaggapadj * repcont
gen union1 = d.union*parcont
gen dunion = d.union*demcont
gen runion = d.union*repcont
gen massinit = massecon*init1
gen mingap = d.newmingap 
gen lmingap = l.mingap
gen avinc = aveincadj * 1000000
gen lavinc = l.avinc
gen nstates = d.minwag2adjneighave2 
gen lnstates = l.nstates
gen init = init1
gen ldemcont = l.demcont
gen lrepcont = l.repcont
gen ldmingap = l.dmingap
gen ldunion = l.dunion
gen ldunemp = l.dunemp
gen ldinflat = l.dinflat
gen lrmingap = l.rmingap
gen lrunion = l.runion
gen lrunemp = l.runemp
gen lrinflat = l.rinflat
gen union12 = union
drop union
gen union = d.union12
drop union12
gen lunion = l.union
gen lunemp = l.unemp
gen nonwhite1 = nonwhite
drop nonwhite
gen nonwhite = d.nonwhite1
drop nonwhite1
gen lnonwhite = l.nonwhite
gen lmassecon = l.massecon
gen lftop1 = l.ftop1
gen parcontm = parcont*mingap
gen parcontu = parcont*unemp
gen parconti1 = parcont*inflat
gen parcontun = parcont*union
gen mingap2 = mingap * mingap
gen dmingap2 = demcont * mingap2
gen union2 = union * union
gen dunion2= demcont * union2
gen inflat2 = inflat * inflat
gen dinflat2 = demcont * inflat2
gen unemp2 = unemp * unemp
gen dunemp2 = demcont * unemp2
gen rmingap2 = repcont * mingap2
gen runion2 = repcont * union2
gen rinflat2 = repcont * inflat2
gen runemp2 = repcont * unemp2
gen mingap3 = mingap2 * mingap
gen dmingap3 = demcont * mingap3
gen union3 = union2 * union
gen dunion3= demcont * union3
gen inflat3 = inflat * inflat2
gen dinflat3 = demcont * inflat3
gen unemp3 = unemp2 * unemp
gen dunemp3 = demcont * unemp3
gen rmingap3 = repcont * mingap3
gen runion3 = repcont * union3
gen rinflat3 = repcont * inflat3
gen runemp3 = repcont * unemp3
gen lparcont = l.parcont
gen lparcontm = l.parcontm
gen lnewmingap=l.newmingap
gen demcontnewgap=demcont*newmingap
gen repcontnewgap=repcont*newmingap 
gen parcontnewgap=parcont*newmingap
gen ldemcontnewgap=l.demcontnewgap
gen lrepcontnewgap=l.repcontnewgap
gen lparcontnewgap=l.parcontnewgap

*Table 1 - 1968-2018
*Minimum Wage - Autoregressive Model - Given the d.v. is in levels (rather than dichotomous)I think we need the lagged d.v. in the equation 
*1968-2018
 xtreg rminwag l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage - 
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - 
xtprobit  rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*Appendix A contains data sources rather than statistical analysis

*Appendix B - Bivariate Analysis

xtreg rminwag l.rminwag parcont if year>1967 & year<2019, fe vce(robust)
xtprobit prevail3 parcont if year>1967 & year<2019, pa vce(robust)
xtprobit  rtw3 parcont if year>1967 & year<2019, pa vce(robust) 

*Appendix C - Two Measures of Legislative Professionalism - MDS1 & MDS2
xtreg rminwag  parcont massecon mds1 union unemp inflat nonwhite init year if year>1967, fe vce(robust)
xtprobit prevail3 parcont massecon mds1 unemp inflat nonwhite union init year if year>1967, pa vce(robust)
xtprobit rtw3 parcont massecon mds1 unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

xtreg rminwag  parcont massecon mds2 union unemp inflat nonwhite init year if year>1967, fe vce(robust)
xtprobit prevail3 parcont massecon mds2 unemp inflat nonwhite union init year if year>1967, pa vce(robust)
xtprobit rtw3 parcont massecon mds2 unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*Appendix D - Influences of State Wage Policy without Initiative States as a Control

xtreg rminwag l.rminwag parcont massecon union unemp inflat nonwhite year if year>1967, fe vce(robust)
xtprobit prevail3 parcont massecon unemp inflat nonwhite union year if year>1967, pa vce(robust)
xtprobit rtw3 parcont massecon unemp inflat nonwhite union year if year>1967, pa vce(robust) 

*Appendix F
*Policy Change Results for the Prevailing Wage and Right-to-Work
*Prevailing Wage - 
xtprobit prevail1 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - 
xtprobit  rtw1 parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*****************************************************************************************

*ALTERNATIVE SPECIFICATIONS
*Lagged Party Control Variables
*Minimum Wage - Autoregressive Model - Given the d.v. is in levels (rather than dichotomous)I think we need the lagged d.v. in the equation 
*1968-2018
 xtreg rminwag l.rminwag l.parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage - 
xtprobit prevail3 l.parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - 
xtprobit  rtw3 l.parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*Prevailing Wage and RTW - Population Average, Random Effects and Fixed Effects With and Without Robust Standard Errors
*Prevailing Wage - logit, pa and robust st. errors
xtlogit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Prevailing Wage - logit, pa and non-robust st. errors
xtlogit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa 
*Prevailing Wage - logit, random effects and robust st. errors
xtlogit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967,  vce(robust)
*Prevailing Wage - logit, random effects and non-robust st. errors
xtlogit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967,  
*Prevailing Wage - probit, pa and robust st. errors
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Prevailing Wage - probit, pa and non-robust st. errors
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa 
*Prevailing Wage - probit, random effects and robust st. errors
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967,  vce(robust)
*Prevailing Wage - probit, random effects and non-robust st. errors
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967,  
*Prevailing Wage - logit, fixed effects and non-robust st. errors - Stata doesn't permit robust standard errors with xtlogit fixed effects 
xtlogit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, fe 
*RTW - logit, pa and robust st. errors
xtlogit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*RTW - logit, pa and non-robust st. errors
xtlogit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa 
*RTW - logit, random effects and robust st. errors
xtlogit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967,  vce(robust)
*RTW - logit, random effects and non-robust st. errors
xtlogit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967,  
*RTW - probit, pa and robust st. errors
xtprobit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*RTW - probit, pa and non-robust st. errors
xtprobit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa 
*RTW - probit, random effects and robust st. errors
xtprobit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967,  vce(robust)
*RTW - probit, random effects and non-robust st. errors
xtprobit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967,  
*RTW - logit, fixed effects and non-robust st. errors - Stata doesn't permit robust standard errors with xtlogit fixed effects 
xtlogit rtw3 parcont massecon unemp inflat nonwhite union init year if year>1967, fe 

*Just Political Variables 
*1968-2018
 xtreg rminwag l.rminwag parcont year if year>1967, fe vce(robust)
*Prevailing Wage - 
xtlogit prevail3 parcont year if year>1967, pa vce(robust)
*Right-to-Work - 
xtlogit  rtw3 parcont year if year>1967, pa vce(robust) 

*Panel Matching
teffects psmatch (minwag2adj) (demcont massecon union unemp inflat nonwhite init year)
teffects psmatch (minwag2adj) (repcont massecon union unemp inflat nonwhite init year)
teffects psmatch (prevail3) (demcont massecon union unemp inflat nonwhite init year)
teffects psmatch (prevail3) (repcont massecon union unemp inflat nonwhite init year)
teffects psmatch (rtw3) (demcont massecon union unemp inflat nonwhite init year)
teffects psmatch (rtw3) (repcont massecon union unemp inflat nonwhite init year)
*Alternative Panel Matching 
*Panel Matching
teffects psmatch (minwag2adj) (demcont massecon union unemp inflat nonwhite init year, probit), atet
teffects psmatch (minwag2adj) (repcont massecon union unemp inflat nonwhite init year, probit), atet
teffects psmatch (prevail3) (demcont massecon union unemp inflat nonwhite init year, probit), atet
teffects psmatch (prevail3) (repcont massecon union unemp inflat nonwhite init year, probit), atet
teffects psmatch (rtw3) (demcont massecon union unemp inflat nonwhite init year, probit), atet
teffects psmatch (rtw3) (repcont massecon union unemp inflat nonwhite init year, probit), atet

*Minimum Wage Equation with just unit (state) fixed effects - thus no year fixed effects
xtreg rminwag lrminwag parcont massecon union unemp inflat nonwhite init if year>1967, fe vce(robust)

*Minimum Wage Equation with just year effects - thus no state fixed effects
xtreg rminwag lrminwag parcont massecon union unemp inflat nonwhite init year if year>1967, vce(robust)

*Minimuum Wage - Random Effects - thus no state or year fixed effects
xtreg rminwag lrminwag parcont massecon union unemp inflat nonwhite if year>1967, vce(robust)

*Prevailing Wage - without year fixed effects - thus no state or year fixed effects - next equation doesn't achieve convergence
*xtprobit prevail3 parcont massecon unemp inflat nonwhite union init if year>1967, pa vce(robust)

*Right-to-Work - without year fixed effects
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init if year>1967, pa vce(robust) 

*Since you can't get state fixed effects in xtprobit switching to xtlogit which allows state and year fixed effects but not either robust standard error or population average (pa) in probit equations which is what was recommended to us - see footnote in paper  
*Prevailing Wage - with state and year fixed effects but not pa or robust standard errors
xtlogit prevail3 parcont massecon ftop1 unemp inflat nonwhite union init year if year>1967,  fe 
*Right-to-Work -  with state and year fixed effects but not pa or robust standard errors - Stata won't estimate rtw version because it can't get convergence
*xtlogit rtw3 parcont massecon ftop1  unemp inflat nonwhite union  init year if year>1967, fe 

*Adding Legislative Professionalism - mds1
*Table 1 - Minimum Wage
*1968-2018
 xtreg rminwag  parcont massecon mds1 union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage -  
xtprobit prevail3 parcont massecon mds1 unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work
xtprobit rtw3 parcont massecon mds1 unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*Adding Legislative Professionalism - mds2
*Table 1 - Minimum Wage
*1968-2018
xtreg rminwag  parcont massecon mds2 union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage -  
xtprobit prevail3 parcont massecon mds2 unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work
xtprobit rtw3 parcont massecon mds2 unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*RESULTS FOR THREE SEPARATE ERAS - 1968-2018, 1968-1980 & 1981-2018
*Table 1 - Minimum Wage
*1968-2018
 xtreg rminwag  parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*1968-1980
 xtreg rminwag  parcont massecon  union unemp inflat nonwhite init year if year>1967 & year<1981, fe vce(robust)
*1981-2018
 xtreg rminwag  parcont massecon union unemp inflat nonwhite init year if year>1967 & year>1980, fe vce(robust)
 *Table 1 - Minimum Wage - Autoregressive Model
*1968-2018
 xtreg rminwag  l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*1968-1980
 xtreg rminwag  l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967 & year<1981, fe vce(robust)
*1981-2018
 xtreg rminwag  l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967 & year>1980, fe vce(robust)
 *Table 2 - Prevailing Wage 
*Column 1 - 1968-2018
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*estimates store mw4
*Column 2 - 1968-1980
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967 & year<1981, pa vce(robust)
*estimates store mw5
*Column 3 - 1981-2018
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1980, pa vce(robust)
*estimates store mw6
*esttab mw4, mw5, mw6, se replace stats(r2 r2_a r2_p aic bic N  ll cmd) starlevels(* .10 ** .05 *** .01) 
*Table 3 - Right-To-Work
*Column 1 - 1968-2018
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 
*estimates store mw7
*Column 2 - 1968-1980
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1967 & year<1981, pa vce(robust) 
*estimates store mw8
*Column 3 - 1981-2018
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1980, pa vce(robust) 
*estimates store mw9
*esttab mw7, mw8, mw9, se replace stats(r2 r2_a r2_p aic bic N  ll cmd) starlevels(* .10 ** .05 *** .01) 

*Regression Discontinuity - Democratic Lower House Strength 45%-55% (upper house any percentage)
keep if lhdemp >.4499 & lhdemp <.55001
*Minimum Wage - Autoregressive Model - Given the d.v. is in levels (rather than dichotomous)I think we need the lagged d.v. in the equation 
*1968-2018
xtreg rminwag l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage - won't run with state fixed effects so just year fixed effects
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - won't run with state fixed effects so just year fixed effects
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 
clear
use E:\PartyPoliticsDataset.dta 
set more off
xtset stnum year, yearly
drop if year<1963
drop if stnum==27
*Interacting Macroeconomic Outcomes and Minimum Wage Gap with Party Control
gen rminwag = d.minwag2adj
gen lrminwag = l.rminwag
gen inflation1 = inflation * parcont
gen inflat = inflation
gen linflat = l.inflation
gen dinflat = inflation * demcont
gen rinflat = inflation * repcont
gen unemp1 = unemp * parcont
gen dunemp = unemp * demcont
gen runemp = unemp * repcont
gen aveincadj1 = aveincadj * parcont
gen aveincadj2 = aveincadj * demcont
gen aveincadj3 = aveincadj * repcont
gen mingap1 = d.newmingap * parcont
gen dmingap = d.minwaggapadj * demcont
gen rmingap = d.minwaggapadj * repcont
gen union1 = d.union*parcont
gen dunion = d.union*demcont
gen runion = d.union*repcont
gen massinit = massecon*init1
gen mingap = d.newmingap 
gen lmingap = l.mingap
gen avinc = aveincadj * 1000000
gen lavinc = l.avinc
gen nstates = d.minwag2adjneighave2 
gen lnstates = l.nstates
gen init = init1
gen ldemcont = l.demcont
gen lrepcont = l.repcont
gen ldmingap = l.dmingap
gen ldunion = l.dunion
gen ldunemp = l.dunemp
gen ldinflat = l.dinflat
gen lrmingap = l.rmingap
gen lrunion = l.runion
gen lrunemp = l.runemp
gen lrinflat = l.rinflat
gen union12 = union
drop union
gen union = d.union12
drop union12
gen lunion = l.union
gen lunemp = l.unemp
gen nonwhite1 = nonwhite
drop nonwhite
gen nonwhite = d.nonwhite1
drop nonwhite1
gen lnonwhite = l.nonwhite
gen lmassecon = l.massecon
gen lftop1 = l.ftop1
gen parcontm = parcont*mingap
gen parcontu = parcont*unemp
gen parconti1 = parcont*inflat
gen parcontun = parcont*union
gen mingap2 = mingap * mingap
gen dmingap2 = demcont * mingap2
gen union2 = union * union
gen dunion2= demcont * union2
gen inflat2 = inflat * inflat
gen dinflat2 = demcont * inflat2
gen unemp2 = unemp * unemp
gen dunemp2 = demcont * unemp2
gen rmingap2 = repcont * mingap2
gen runion2 = repcont * union2
gen rinflat2 = repcont * inflat2
gen runemp2 = repcont * unemp2
gen mingap3 = mingap2 * mingap
gen dmingap3 = demcont * mingap3
gen union3 = union2 * union
gen dunion3= demcont * union3
gen inflat3 = inflat * inflat2
gen dinflat3 = demcont * inflat3
gen unemp3 = unemp2 * unemp
gen dunemp3 = demcont * unemp3
gen rmingap3 = repcont * mingap3
gen runion3 = repcont * union3
gen rinflat3 = repcont * inflat3
gen runemp3 = repcont * unemp3
gen lparcont = l.parcont
gen lparcontm = l.parcontm
gen lnewmingap=l.newmingap
gen demcontnewgap=demcont*newmingap
gen repcontnewgap=repcont*newmingap 
gen parcontnewgap=parcont*newmingap
gen ldemcontnewgap=l.demcontnewgap
gen lrepcontnewgap=l.repcontnewgap
gen lparcontnewgap=l.parcontnewgap
keep if uhdemp >.4499 & uhdemp <.55001
*Democratic Upper House Strength Between 45%-55% (lower house any percentage)
keep if lhdemp >.4499 & lhdemp <.55001
*Minimum Wage - Autoregressive Model - Given the d.v. is in levels (rather than dichotomous)I think we need the lagged d.v. in the equation 
*1968-2018
xtreg rminwag l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage - won't run with state fixed effects so just year fixed effects
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - won't run with state fixed effects so just year fixed effects
xtprobit rtw3 parcont massecon unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 
*Democratic Strength in Both Champers Between 45%-55%
keep if lhdemp >.4499 & lhdemp <.55001
*Minimum Wage - Autoregressive Model - Given the d.v. is in levels (rather than dichotomous)I think we need the lagged d.v. in the equation 
xtreg rminwag l.rminwag parcont massecon union unemp inflat nonwhite init year if year>1967, fe vce(robust)
*Prevailing Wage - won't run with state fixed effects so just year fixed effects
xtprobit prevail3 parcont massecon unemp inflat nonwhite union init year if year>1967, pa vce(robust)
*Right-to-Work - won't run with state fixed effects so just year fixed effects
xtprobit rtw3 parcont massecon  unemp inflat nonwhite union  init year if year>1967, pa vce(robust) 

*************************************************************************************************************************************
*Appendix E - See discussion at the beginning of this do file.

*use E:\PartyPoliticsDataset.dta
*set more off
*drop id
*gen id=stnum
*encode id, gen(nid)
*xtset id year, yearly
*clustergen rminwag prevail3 rtw3 parcont massecon union unemp inflat nonwhite init year 

*gen lbwrminwag_bw = l.rminwag_bw 
*gen lwirminwag_wi = l.rminwag_wi
*gen lbwprevail3_bw = l.prevail3_bw  
*gen lwiprevail3_wi = l.prevail3_wi
*gen lbwrtw3_bw = l.rtw3_bw  
*gen lwirtw3_wi = l.rtw3_wi
*gen lbwparcont_bw = l.parcont_bw  
*gen lwiparcont_wi = l.parcont_wi
*gen lbwmassecon_bw = l.massecon_bw  
*gen lwimassecon_wi = l.massecon_wi
*gen lbwunion_bw = l.union_bw  
*gen lwiunion_wi = l.union_wi
*gen lbwunemp_bw = l.unemp_bw  
*gen lwiunemp_wi = l.unemp_wi
*gen lbwinflat_bw = l.inflat_bw  
*gen lwiinflat_wi = l.inflat_wi
*gen lbwnonwhite_bw = l.nonwhite_bw  
*gen lwinonwhite_wi = l.nonwhite_wi
*gen lbwinit_bw = l.init_bw  
*gen lwiinit_wi = l.init_wi
*gen lbwyear_bw = l.year_bw  
*gen lwiyear_wi = l.year_wi

*xtreg  rminwag l.rminwag parcont_wi massecon_wi union_wi unemp_wi inflat_wi nonwhite_wi init_wi parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi     if year>1967, i(id) mle
*xtreg  rminwag l.rminwag parcont massecon union unemp inflat nonwhite init parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi if year>1967, i(id) mle
*xtprobit prevail3 l.prevail3 parcont_wi massecon_wi union_wi unemp_wi inflat_wi nonwhite_wi init_wi parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi  if year>1967, pa vce(robust)
*xtprobit prevail3 l.prevail3 parcont massecon union unemp inflat nonwhite init parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi if year>1967, pa vce(robust)
*xtprobit rtw3 l.rtw3 parcont_wi massecon_wi union_wi unemp_wi inflat_wi nonwhite_wi init_wi parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi  if year>1967
     
*xtprobit rtw3 l.rtw3 parcont massecon union unemp inflat nonwhite init parcont_bw  massecon_bw  union_bw  unemp_bw  inflat_bw  nonwhite_bw init_wi if year>1967
